A Re normalization Group computation of the critical exponents of hierarchical 

spin glasses 



Michele Castellana^'^ and Giorgio Parisi^ 

^ Dipartimento di Fisica, Universita di Roma 'La Sapienza' , 00185 Rome, Raly 
^LPTMS, CNRS and Umversite Pams-Sud, UMR8626, Bat. 100, 91405 Orsay, France 

In a recent work (M Castellana and G Parisi, Phys. Rev. E 82, 040105(R) (2010)), the large- 
scale behaviour of the simplest non-mean field spin-glass system has been analysed, and the critical 
exponent related to the divergence of the correlation length computed at two loops within the 
e-expansion technique with two independent methods. 

By performing the explicit calculation of the critical exponents at two loops, one obtains that the 
two methods yield the same result. This shows that such underlying renormalization group ideas 
apply also in this disordered model, in such a way that an e-expansion can be consistently set up. 
The question of the extension to high-orders of this e-expansion is particularly interesting from the 
physical point of view. Indeed, once high orders of the series in e for the critical exponents are 
known, one could check the convergence properties of the series, and find out if the ordinary series 
re-summation techniques yielding very accurate predictions for the Ising model work also for this 
model. If this is the case, a consistent and predictive non-mean field theory for such disordered 
system could be established. 

In that regard, in this work we expose the underlying techniques of such a two-loop computation 
(M Castellana and G Parisi, Phys. Rev. E 82, 040105(R) (2010)). We show with an explicit 
example that such a computation could be quiet easily automatized, i. e. performed by a computer 
program, in order to compute the e-expansion at high orders in e, and so eventually make this theory 
physically predictive. Moreover, all the underlying renormalization group ideas implemented in such 
a computation are widely discussed and exposed. 
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I. INTRODUCTION 



The understanding of glassy systems and their 
critical properties is a subject of main interest in 
Statistical Physics. The mean-field theory of spin- 
glasses [ij and structural glasses provides a phys- 
ically and mathematically rich theory. Neverthe- 
less, real Spin-Glass systems have short-range inter- 
actions, and thus cannot be successfully described 
by mean- field models This is the reason why 
the development of a predictive and consistent the- 
ory of glassy phenomena going beyond mean field 
is still one of the most hotly debated, difficult and 
challenging problems in this domain 0-0], so that 
a theory describing real glassy systems is still miss- 
ing. Indeed, the standard field theory techniques 
0jQ yielding the Ising model critical exponents with 
striking agreement with experimental data do not 
usually apply to locally-interacting glassy systems. 
As a matter of fact, a considerable difficulty in the 
set-up of a loop-expansion for a locally-interacting 
Spin-Glass is that the mean-field saddle-point has 
a very complicated structure 0], and could be not 
uniquely-defined It follows that the predictions 
of a loop-expansion performed around one selected 
saddle-point could actually depend on the choice of 
the saddle-point itself, resulting into an intrinsic am- 



biguity in the physical predictions. Moreover, non- 
perturbativc effects arc poorly understood and not 
under control, and the basic properties of large scale 
behaviour of these systems arc still far from being 
clarified. 

In ferromagnetic systems, the physical properties 
of the paramagnetic- ferromagnetic transition emerge 
in a clear way already in the original approach of 
Wilson Q, where one can write a simple Renormal- 
ization Group (RG) transformation. It was later 
realized that Wilson's equations are exact in the in 
models with ferromagnetic power-law interactions 
on hierarchical lattices as the Dyson model [1, • 
Indeed, in such models one can write exact equa- 
tions for the magnetization probability distribution, 
containing all the relevant physical informations 
about the paramagnetic, ferromagnetic and critical 
fixed point, and the existence of a finite-temperature 
phase transition. In other words, all the physical 
RG ideas are encoded in such recursion relations, 
whose solution can be explicitly built up with the 
e-expansion technique [T]| . 

The extension of this approach to random systems 
is available only in a few cases. 

Firstly, an RG analysis for random models on the 
Dyson hierarchical lattice has been pursued in the 



past [12, ll3| , and a systematic analysis of the phys- 
ical and unphysical infrared (IR) fixed points has 
been developed within the e-expansion technique. 
Unfortunately, in such models spins belonging to the 
same hierarchical block interact each other with the 
same [l^l random coupling J, in such a way that 
frustration turns out to be relatively weak and they 
arc not a good representative of realistic strongly 
frustrated system. 

Secondly, models with local interactions on hier- 
archical lattices built on diamond plaquettes [23j . 
have been widely studied [23 - l28j in their spin glass 
version, and lead also to weakly frustrated systems 
even in their mean- field limit [29[. Notwithstanding 
this, such models yield a very useful and interesting 
playground to show how to implement the RG ideas 
in disordered hierarchical lattices, and in particular 
on the construction of a suitable decimation rule for 
a frustrated system, which is one of the basic topics 
in the construction of a RG for spin-glasses, and so 
in the identification of the existence of a spin-glass 
phase in finite dimension. 

Moreover, there has recently been a new wave of 
interest for strongly frustrated random models on 
hierarchical lattices (T^ - [l6j : for example, it has been 
shown that a generalization of the Dyson model 
to its disordered version (the Hierarchical Random 
Energy Model (HREM)) has a Random Energy 
Model-like phase transition, yielding interesting 
new critical properties that don't appear in the 
mean-field case. 

In a recent work [30j . we performed a field theory 
analysis of the critical behavior of a generalization of 
Dyson's model to the disordered case, known as the 
Hierarchical Edwards- Anderson model (HEA) 
that is physically more realistic than the HREM 
and presents a strongly-frustrated non-mean field 
interaction structure, being thus a good candidate 
to mimic the critical properties of a real spin-glass. 
Moreover, the symmetry properties of the HEA 
make an RG analysis simple enough to be per- 
formed with two independent methods, to check if 
the IR-limit of the model is physically well-defined 
independently on the computation technique that 
one uses. Another element of novelty of the HEA 
is that its hierarchical structure makes the RG 
equations simple enough to make an high-order 
e-expansion tractable by means of a symbolic 
manipulation program, resulting in a quantitative 
theory for the critical exponents beyond mean field 
for a strongly-frustrated spin-glass system. It is 
possible that such a perturbative expansion turns 
out to be non convergent: if this happens, it may 



help us to pin down the non-perturbative effects. 
Motivated by this purpose, we have shown js^l with 
a two-loop calculation that such e-cxpansion can 
be set up consistently, and that the ordinary RG 
underlying ideas actually apply also in this case, 
so that the IR limit of the theory is well-defined 
independently on the regularization technique. 

In the present work, we show how the underlying 
RG ideas emerge in the computation, and in partic- 
ular how such a calculation has been performed, so 
that the reader can fully understand and reproduce 
it. Moreover, we show by an explicit example of 
such computation how this e-expansion could be 
automatized, i. e. implemented by a computer 
program, and so pushed to high orders to establish 
its summability properties. 

The HEA is defined [H, (M] as a system of 2'^+^ 
Ising spins Sq, . . . , S'2fc+i_i, with an energy function 
defined recursively by coupling two systems of 2*^ 
Ising spins 

^fc+i [•5'o, ■ • ■ , 5'2fc+i-i] = (1) 
H'l^ [So, S'2fc_i] + H'l^ [S2k, . . . , 52/.-+i_i] + 



0,2 



2{k+l)a 



J ij S[ Sj , 



where 



i?/[S'i, 5*2] — — J2 '^8182 



and Jij, J are Gaussian random variables with zero 
mean and unit variance. 

As we will show in the following, this form of the 
Hamiltonian corresponds in dividing the system in 
hierarchical embedded blocks of size 2*^ and that the 
interaction of two spins depends on the distances 
of the blocks to which they belong [l3, HI]- Here 
< cr < 1 is a parameter tuning the decay of 
the interaction strength with distance. The HEA 
is therefore a hierarchical counterpart of the one- 
dimensional spin glass with power-law interactions 
which has received attention recently [32|-|36|. 

It is crucial to observe [lj| that the sum of the 
squares of the interaction terms that couple the two 
subsystems scales with k as 2^'^(^~°'). Hence, for 
(T > 1/2 the interaction energy scales sub-extensively 
in the system volume, yielding a non-mean field be- 
havior of the model, while for a < 1/2 it grows faster 
than the volume, and the thermodynamic limit is not 
defined. The interesting region we will study is thus 
a e [1/2,1]. 

An equivalent definition of the HEA can be given 



without using the recursion relation ([T]). Indeed, one 
can recover ((T]) by defining the HEA as a system of 
2^ Ising spins with Hamihonian 

where J^- are Gaussian random variables with zero 
mean and variance a^y The form of af^ is given 
by the following expression: if only the last m 
digits in the binary representation of the points i 
and i are different, af,^ — 2~^'^™. This form of the 
Hamiltonian corresponds in dividing the system in 
hierarchical embedded blocks of size 2'", such that 
the interaction between two spins depends on the 
distance of the blocks to which they belong. It is 
important to observe that the quantity cr^^- is not 
translational invariant, but it is invariant under a 
huge symmetry group and this will be crucial in the 
study of the model. 

The two definitions ([1]) and ([2]) are equivalent. 



We reproduce the IR behavior of the HEA and 
calculate its critical exponent by two different meth- 
ods. Both methods suppose the existence of a grow- 
ing correlation scale length 5, diverging for T — > Tc 
as 

^(x(T-T,)-^ 

in such a way that for T — > Tc the theory is invariant 
under rc-parametrizations of the length scale. 
The first method is analogous to the coarse-graining 
Wilson's method for the Ising model: the scale- 
invariant limit is obtained by imposing invariance 
with respect to the composition operation of Eq. 
([1]), taking two systems of 2'^ spins and yielding a 
system of 2*^+^ spins. As for the Dyson ferromag- 
netic model, thanks to the hierarchical structure of 
the Hamiltonian one can obtain closed formulae for 
physical quantities with respect to such composition 
operation, analyse the critical and non-critical 
fixed-points and extract v. 



The second method is more conventional: the IR 
divergences appearing for T ^ Tc are removed by 
constructing a renormalized IR-safe theory. The 
fundamental physical informations one extracts 
from such renormalized theory are the same as 
those of the original theory defined by Eq. ([1]). In 
particular, the correlation length and its power-law 
behavior close to the critical point must be the 
same, and so the critical exponent ly. 

The rest of this paper is divided into three main 
Sections: in Sec. |II]we go through the main steps of 
the computation with Wilson's method, show that 
the tensorial operations can be easily implemented 
diagrammatically, and thus performed by a com- 
puter program to extend such e-expansion to high 
orders. Moreover, we give the two-loop result for i>. 
In Sec. mil the same result is reproduced with the 
field-theoretical method, and the analogies between 
the two methods are discussed. In particular, we dis- 
cuss why Wilson's method would be definitely better 
for an automatization of the e-expansion to high or- 
ders. Both in Sections |TT] and IIIIl we explicitly do 
all the steps of the calculation at one loop, giving to 
the reader all the information needed to reproduce 
the two- loop result for v. 

Finally, in Sec. llVl the two- loop result is discussed 
in the perspective of the set up of a high-order e- 
expansion. 



II. WILSON'S METHOD 

As mentioned before, the hierarchical symmetry 
structure of the model makes the implementation of 
a recursion-like RG equation simple enough to be 
solved within an approximation scheme. As a matter 
of fact, let us define the probability distribution of 
the overlap P, 

Qab, a = 1, ■ • • ,n, (3) 

Qab = Qba, Qaa = Vfl, = 1 , . . . , n 

as 



2"-l 



= Ei 



E exp -/3 E Hi [^o", ■ • • , ^2^-1] n M ^"fc - 4 E 



(4) 



\{S^)^ 



a=l 



a<b=l 



and the rescaled overlap-distributions as 



where /3 = 1/T is the inverse temperature. 
According to the general prescriptions of the replica 



approach, all the Physics of the model is encoded in 
the n -> limit of Zk[Q]. 

It is easy to show that the recursion relation ((!]) 
for the Hamiltonian results into a recursion relation 
forZfe[Q] 



Zfc[Q]=exp(^Tr[Q2] 



(5) 



[dP] Z, 



k-l 



Q + P 



X Z, 



Q-P 



where Tr denotes the trace over the replica indexes, 
/ [dP] stands for the functional integral over Pab and 



C = 2^(1- 



(6) 



Eq. ([S]) physically represents a recursion equation 
relating the probability distribution Zk-i[Q] to 
Zk[Q], obtained by Zi;~i[Q] by a coarse-graining 
RG step, composing two subsystem of size 2*^ to 
form a system of size 2'^+^. Eq. ([5]) is analo gou s 
to the recursion equation in Dyson's model 
relating the probability distribution 5fc(m) of the 
magnetization at the fc-th hierarchical level to 
gk-i{m). 



To illustrate the technique used to solve pertur- 
batively (0) for Zk[Q], we will show our method in 



a simple toy example, where the matrix field Qab is 
replaced by a one-component field 0, the functional 
Zk[Q] by a function rife (</>), and Eq. ^ by 



Qki4>) = cxp 
dx^k-i 



Cl/2 



(7) 



X O 



fc-i 



(f>-x 

Cl/2 



As for Dyson's model, ([7]) can be solved by mak- 
ing an ansatz for fife (</)). The simplest form one can 
suppose for rife(</') is the Gaussian one 



17fe(0) =exp [-(40^ + efe(^)] 



(8) 



This form corresponds to a mean- field solution 
By inserting Eq. ^ into Eq. d?]), one finds a recur- 
sion equation relating 4, efe to 4-1, efe-i 



dk = 



2d, 



fc-i 



C 
2efc-i 
Cl/2 



4 ' 



Non-gaussian solutions can be explicitly constructed 
perturbatively. Indeed, by setting 



Zk[Q] = exp - ( dkc/) + ek4> + -5- 



Uk ,3 



(9) 



and supposing that Uk is small, one can plug Eq. ([9]) 
into Eq. ([7]) and get 



VLk{(t)) = exp 

oc exp 

= cxp 

oc exp 

oc exp 
1 

+ 3 




(10) 



r 



where oc stands for a 0-indepcndcnt proportionality 
constant, that will be omitted in the following. Com- 
paring Eq. (fin)) with Eq. ©, one finds three recur- 



rence equations relating 4, &k, Uk to 4-i, efc-i, Uk-i 
,2 . / N 2 



24-1 



C 



Uk-1 



4 4 \2C^/^dk-i 



(11) 



2e 



efc = 



Uk 



Ci/2 2Ci/24._i 
C3/2 + 2 I 2Ci/2d 



One can easily analyse the fixed points of the 
RG-flow equations pT|) . and the resulting critical 
properties. We will not enter into these details for 
the toy model, since all these calculations will be 
illustrated extensively in the original theory. 



Back to the original problem, Eq. ^ can be 
solved by making an ansatz for Zk [Q] , following the 
same lines as in the toy model case. The simplest 
form one can suppose for Zk [Q] is the Gaussian one 



cxp (-rfeTr[Q2]) 



(12) 



This form corresponds to a mean- field solution Q • 
By inserting Eq. (fT2|) into Eq. (O, one finds the 
evolution equation relating rk-i to rj. 



rk 



C 



4 ' 



(13) 



Corrections to the mean-field solution can be inves- 
tigated by adding non-Gaussian terms to Eq. (fT2|) . 
that are proportional to higher powers of Q, and con- 
sistent with the symmetry properties of the model. 
It is easy to see [l| that the only cubic term in Q 
consistent with such symmetry conditions is Tr[(3'^], 
so that the non-mean field ansatz of Z reads 



cxp 



Wk 



(14) 



This correction can be handled by supposing that 
the term Wk proportional to the non-quadratic term 
is small for every k, and performing a systematic 
expansion in powers of it. By inserting Eq. p4p into 
the recursion relation Eq. ((51), one finds 



^"^-iTr [P2] + ^!^Tr[gp2i 



C 



C3/2 



2wk- 



3C3/2 



iTr[0^ 



[dP] exp 



(15) 



The Gaussian integral in Eq. p3|) can be computed Mab,cd[Q] = Nab,cd[Q] + Nab,dc[Q], (17) 

exactly. Indeed, setting A = {a,b)a>b and N^hMQ] = ^bcQda + SacQdb, (18) 

—^^ab + ^^^^^Mab[Q], one finds 



Zk[Q] = cxp <^ - 



2rk- 
C 



4 



Tr [Q2 



2wfc- 
3C3/ 



iTr[Q3 



det 



4r 



k-l 



C 



JAB 



-Mab[Q] 



(19) 



The determinant in the right-hand side of Eq. (fTSj) evaluate the traces Tr[M[Q]^], Tr[Af [Q]^]. Here we 
can now be expanded in Wk-i- Denoting by Tr the show how the trace Tr[M[Q]^] can be evaluated, in 
trace over the A-type indexes, one has to explicitly order to show to the reader how the tensorial opera- 



tions over the replica indexes can be generally carried By using Eqs. (|17p . ([T5)l . one has 
out. 



I 

Tr[M[Qf] = J2m[Q]abM[Q]ba (20) 

AB 

= J2 (NabMQ] + NabMQ]) X 
a>b.c>d 

X iNcd,ab[Q] + Ndc^abiQ]) 
= J2 Nab,cd[Q]Ncd,ab[Q] 

= X] i^bcQda + 5acQdb) X 
a^b^c^d 

X {SdaQbc + ^caQbd) 

= X ^caQbd 
a^b^c^d 

= ^(l-M(l-M^caQL 
abed 

= (n-2)Tr[Q2]. 

I 



The steps in Eq. (|20p can be summarized as follows: 
in the second line we write the tensorial operations 
over the super indexes A,B, -- in terms of the 
replica indexes a,b, - ■ ■ , in the third line we use the 
symmetries of Mab.cd[Q] with respect to a o 6 and 
c d and re-write the sum over a > b, c > d in 
terms of a sum with a b, c ^ d. In the fifth line we 
find out which amongst the terms stemming from 
the product {SbcQda + SacQdb)iSdaQbc + ScaQbd) van- 
ish because of the constraints a ^ b, c ^ d, Qaa ~ 0, 
and because of the Kronecker 5s in the sum. Once 
we are left with the non- vanishing terms, in the sixth 
line we write explicitly the sum over a ^ b, c ^ d 
in terms of an unconstrained sum over a, 6, c, d by 
adding the constraints (1 — 5ab)(l — Scd)- In the 
seventh line we perform explicitly the sum over the 
replica indexes, and write everything in terms of the 
invariant 4^\Q] = Tt[Q^] (see Table 



Here we show how the trace in Eq. (PD|) can be 
computed with a purely graphical method, that can 
be easily implemented in a computer program to per- 
form this computation at high orders in Wk- Let us 
set 

Tr2[/] EE fa.bu-.akbk (21) 

H (l-^aibi)---(l-'5a,!,J X 
aihi,--- jQfcbfc 

x/aibi,--- .Ofcbfc ■ 

and make the graphical identifications shown in Fig. 

m 



The last line in Eq. ([20]) can be reproduced by 
a purely graphical computation, as shown in Fig. 
[21 There we show that all the tensorial operations 
have a graphical interpretation, and so that they 



can be performed without using the cumbersome 
notation of Eq. pOj) . This graphical notation 
is suitable for an implementation in a computer 
program, that could push our calculation to very 



FIG. 1: Graphical identifications representing symbolically mathematical objects used in tensorial operations. The 
basic objects of the tensorial operations are the Sat function imposing that the replica indexes a and b are equal (top), 
and the matrix Qat (bottom). Once these elements are represented graphically, all the tensorial operations can be 
carried out by manipulating graphical objects composed by these elementary objects. 



high orders. For example, as shown in Fig. [2] in 
a simple example, while computing Tr[Af [Q]*^] for 
fc ^ 1, a proliferation of terms occurs, and some of 
these terms can be shown to be equal to each other, 
because represented by isomorph graphs, so that 
the calculations can be extremely simplified. For 
example, in a computer program implementation of 
such a computation, one could identify the isomorph 
graphs by using some powerful graph isomorphism 
algorithm [3lj |. 



Comparing Eq. ([23]) with Eq. (fT4)) . one finds a re- 
cursion relation for the coefficients , Wk 



Tk = 
Wk = 



Eq. dini) shows that for e = cr - 2/3 < Wfe ^ for 
k oo, i. e. the corrections to the mean-field Zk[Q] 
vanish in the IR limit. In this case, the critical fixed 
point (r*, w*) of Eqs. (|42]), (gO]) has 0. On the 

contrary, for e > a non-trivial critical fixed-point 
7^ arises. According to general RG arguments, 
such non-trivial fixed point will be proportional to e 
0. In particular, one finds that w1 — 0{e). 



By following these steps shown Eq. ([20]) (or their 
graphical implementation) , all the other tensorial op- 
erations can be carried out. In particular, one finds 

Tr[Af[g]3] = (n-2)Tr[03]. (22) 



By plugging Eqs. ^ into Eq. one finds 



(23) 



r 



(24) 



I 

The critical exponent v can be computed Q by 
considering the 2x2 matrix Ai linearising the trans- 
formation given by Eq. ([24]) around the critical fixed- 
point (r*,?!!*) 



ywk-w^ J y Wk^i -w^ J ' 



Zk[Q] = exp< - 



2rk- 



C 



4 



n-2 



Wk-l 



Tr[g2 



1 f 2wk-i 
3 I C3/2 



Wk-l 

2C1/2 



TT[Q'] + Oiwt_,) 



2rk-i 
C 

2wk-i 

C3/2 



4 



n-2 



Wk-l 



+ Oiwt,), 



IC^I^Tk-l 

+ 0(<-_i). 



2C1/2 



rk-i 



and is given by 



where A is the largest eigenvalue of M . 



log 2 
log A' 



(25) 



Such procedure can be systematically pushed to 
higher order in w^, and thus in e, by taking into ac- 
count further corrections to the mean-field solution. 



Indeed, if we go back to Eq. (|19p and consider also 
the 0{w'^_^) terms in the right side, we find 



J 



exp ■ 



-Tr 



2Ci/2r,._i 



C 



C3/2 
■Wk-l 



Wk-1 



4 [ 2C7i/2r,._i 
I 



1 

4 

M[Q]^ + 0{wl^, 



(26) 



By computing explicitly the 0(z/;^_]^) term in the 
right side of Eq. ((26l) , one finds 



Tr[M[QY 



(4)r 





-16/f [Q]- 




= Tr[Q4], 




^ (Tr[Q2])2^ 




— ^ ^ QabQacJ 








^ T.Q'b- 




ab 


Zk[Q] = exp 



(27) 



Plugging Eq. ^ in Eq. and Eq. ^ in 

Eq. (HH), we see that at 0{wl^-^), Eq. © gen- 
erates the fourth-order monomials {I^^^ [Q]}i=i^...a, 
that are not included into the original ansatz ([T4|). 
It follows that at 0{w'^), Zk[Q] must be of the form 



r.Tr[Q^]+^Tr[Q3] + iEA,,/(^)[Q] 



(=1 



(28) 



withA,fc = 0(u;4)V/ = l,...,4. 



By inserting Eq. p8)) into Eq. ^ and 
expanding up to 0{w'j._^), we obtain six re- 
cursion equations relating rk,Wk, Xik, ■ ■ ■ to 
rk-i,Wk-i, Xik-i, ■ ■ ■ ,X4k-i- 

Such systematic expansion can be iterated to any 



order O(w^), obtaining 



p "j 



Z,[Q]=cxp(-E5]c|;l//^")[Q]|, (29) 



where cfl = rk, c\^l = Wfe/3,cJ'^^ = Xikf-i^l = 
1, ... ,4. In this way, a recursion equation relating 



j=2 1=1 
„(3) - ,„ /Q .(4) 



Example of a graphical computation 



Tr[M[Qf] = Tr2[N[Q]N[Q]] 





+ 



2x 




+ 




o 



= Tr 





- 2x 



nx 



= (n - 2)Tr[Q2] 
FIG. 2: Graphical computation of Tr[Af[Q]^] in Eq. (|20|| . 

In the second line, the two addends of the matrix N[Q]ab,cd in Eq. (|18p are represented graphically in terms of the 
graphical objects defined in Fig. [1] 

In the third line, the legs of such addends are contracted with each other, and four terms are generated. Two of them 
are topologically identical, hence we simply put a factor 2 multiplying the second term in the second line. 
According to the condition Qaa = in Eq. (|3}, the first and second term in the third line vanish. Indeed, in these 
terms the lines coming out of the square vertex (Qab) are connected by a circuit, meaning that the matrix element 
Qab is computed with a — b, and thus vanishes. Moreover, the two top-lines in the third term are actually equivalent 
to just one line, because of the relation S^i, = Sab- Hence, we are left with a single term in the fourth line. 
In the fifth line, we perform graphically the operation Tr2. Such an operation can be easily implemented graphically 
by looking at the second line of Eq. (|21|l . Let us expand the product of Ss in the second line of Eq. (I21|) . and recall from 
Fig. [T]that Sab represents a line with a circular dot connecting a with b. Hence, given graphical object O, with external- 
legs (indexes) (ai, ■ ■ ■ , (a^, fe^), Tr2[0] is nothing but a sum of all the possible 2* contractions (performed with a 
line with a circular dot) of such external legs, where each contracted term is muhiplied by (-1)* °' °f th<= 

In this case fc = 2, so we generate 2^ terms in the fifth line. 

In the sixth line, we take into account the fact that the second and third term in the fifth line are topologically 
isomorph, and that the fourth term in the fifth line vanishes because of the relation Qaa = 0. 

In the seventh line the unconstrained sum over the replica indexes is finally performed. This can be done graphically 
in the following way: when we have an external line with a round vertex, summing over the replica index represented 
by that line means that one has simply to remove the line (this is the graphical implementation of the relation 
<5ai)<7i) = Qa)- We do this in the first term: we sum over the top-left index, and remove the line on the top. Then 
we sum over the top right index by simply multiplying by n. The sum over the bottom-left and bottom-right indexes 
simply yields TrlQ'^]. We to the same for the second term: we sum over the top-right index and remove the top line, 
then sum over the top-left index and remove the top-left line. Then, the sum over the bottom-left and bottom-right 
indexes yields Tr[Q^]. 

Hence, we get the same result as in Eq. (|20|l . 



{c; fc-ili,/ to {cjffc}j,i is obtained. 

The number rij of monomials generated at the step 
j of this procedure prohferates for increasing j. In 
Table |I] we show the invariants [Q] obtained by 
performing this systematic expansion up to the or- 
der p = 5. It is interesting to observe that the in- 



variants Tr[Q2]2^ Tr[Q2]Tr[Q3] that are generated, 
are of 0{n^) if the matrix Qab is replica symmet- 
ric. Notwithstanding this, in general they will give a 
non- vanishing contribution to the recursion relations 
Wi fc-ilj.i ~^ Wi k\ 3,1-1 8'iid so to V. The recurrence 
equations at this order are the following 
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-Ik 
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-Ik 
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-Ik 
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^3k 



^4k 
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-Ik 
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^3k 



^4k 
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'-Ik-l 
J4) 
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+0((c(l,)«). 



(4) o (4) 



[12 + n[n - 1)] 



Q„(5) Opi^^^ 
20C3/2cffc^_l 20C3/2cffcLi 
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-3/C-I 
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,(3) 
1 /c-1 



fe-i 



C2 



2c 



.(4) 

-2 fe-1 



2 V2CV242)_^^ 
3 



C2 




+ 0((cll,)«). 



+ 0((cfL,f), 
+ 0((cfL)^ 



((41:)^ 



2c^") 



9/^„(2) OcJ-'^'l 
'^^'^lk-1 ^^^^Ife-l. 



(30) 



(31) 



(32) 
(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 



By looking at Eqs. ((32|) - ([39)) and using the cj^ , C;\, scale to zero as A: -> 00 if e < 1/12. It 
definition it is easy to see that the coefficients 



j 




2 


Tr[Q^] 


3 




4 


Tr[Q*] 


Tr[Q^]^ 


X^aT^c QabQbc 




5 


Tr[Q''] 


Tr[Q"]Tr[Q^] 


J2abcd QabQbcQbdQcd 


X^atc QabQacQbc 



TABLE I: Invariants generated to the order p — 5. In each line of the table we show the invariants [Q] , . . . , I„- [Q] 
from left to right. 



is easy to find out that this is actually true for 
all the coefficients cp^ with j > 3. Such a critical 
value e = 1/12 will be reproduced also in the 
field-thcorctical approach in Section Hill 

The evolution Eqs. ([30]) - (fSQ]) depend smoothly on 
the replica number n, so that the analytical continu- 
ation ri — > 0, can be done directly. By linearising the 
transformation (pO|l - (|39)) around the critical fixed- 
point {cp^lj,; and computing the matrix Ai, one can 
extract A, and so v for 7i = to the order by using 
Eq. (ESI). We find 

f = 3 + 36e -I- [432 - 27(50 + 55-21/3+ (40) 
+53-22/3) log 2] £2 + 0(63). 



The one-loop result for v is the same as that 
of the power-law interaction spin-glass [1] (where 
e = 3((T — 2/3)). Notwithstanding this, the coeffi- 



III. FIELD-THEORETICAL METHOD 

Here the IR limit is performed by constructing 
a functional integral field theory and by removing 
its IR divergences within the minimal subtraction 
scheme. 

While in Wilson's method the IR limit was 
reached by looking at the scale invariant fixed- 
points of the recursion relation ([5]) for fc — > oo, in 
this case the wc will take before the large-fc limit, 
remove the resulting IR singularities through re- 
normalization, and then perform the scale-invariant 
limit by means of the Callan-Symanzik equation. 

This computation is better performed by slightly 
changing the definition of the model. Indeed, the 
following re-definition of the interaction term in Eq. 



cients of the expansion in these two models will be 
in general different at two or more loops. As a matter 
of fact, the binary tree-structure of the interaction of 
the HEA emerges in the non-trivial log 2, 2^/^ factors 
in the coefficient of in (|40|) . that come from the 
binary structure of the hierarchical tree and can't be 
there in the power-law case. 

Before discussing the result in Eq. ((40|) . we point 
out that Wilson's method explicitly implements the 
binary-tree structure of the model when approach- 
ing the IR limit. As a matter of fact, the hierarchi- 
cal structure of the model is explicitly exploited to 
construct the steps of the RG transformation. Nev- 
ertheless, if the IR limit is unique and well-defined, 
physical observables like v must not depend on the 
technique we use to compute them in such a limit. 
It is thus important to verify that Eq. (|40|) does 
not depend on the method we used to reproduce the 
IR behavior of the theory. This has been done by 
reproducing Eq. (|40|) with a quite different field- 
theoretical approach. 

I 

m 

J2 JrjS.S, ^ E E J^^^^^^- (41) 

i<j i=0 j=2'' 

is equivalent to the original definition ([1]) and 
makes the field theory computations simpler. The 
equivalence of (14 It with the original definition ([T]) 
can be shown |1^ by observing that the scaling 
of the spin coupling in the model defined by Eq. 
((4T|) differs from that in Eq. ([1]) for a constant 
multiplicative factor, and thus that the two options 
are equivalent, and must yield the same critical 
exponents. Notwithstanding this, the critical tem- 
perature of the model defined by Eq. ([T|) and that 
of the model defined by Eq. (|4T|) will be different, 
since the spin couplings Jtj differ by a multiplicative 
factor, resulting in a redefinition of /?. Physically 
speaking, the original definition ([1]) is such that. 



when two subsystems of 2^ spins arc coupled to 
form a system with l}^^^ spins, one introduces 
couphngs Jij between the two subsystems and 
additional couplings between the spins within each 
subsystem, while in (j4ip only couplings between the 
two subsystems are introduced. 

By iterating the recursion relation for the Hamilto- 
nian, one has an explicit form for \S\ of a system 
of 2^ spins in the large-A: limit. Then, the average of 
the replicated partition function is be expressed as 
an integral over the local overlap field Qiah = S°;S\ 



Ej 



^ exp -/?^i^,^[5o^...,52"._l] 

{S'°}.,a V a=\ / 



\dQ\ 



By using dimensional analysis, it is easy to pick up 
the terms in ^[Q] that are relevant in the IR-limit. 
It is easy to check that S\Q\ is given by the sum of 
a quadratic term in Qiah-, plus a cubic term, plus 
higher-degree terms. The dimensions of the field 
Qi ah can be computed by imposing the adimension- 
ality of the quadratic term, and so the dimension of 
the coefficient g of the cubic term and of those of 
the higher-degree terms. One finds that the dimen- 
sions of g in energy is \g\ = 3e. Thus, as in Wilson's 
method, the cubic term scales to zero in the IR limit 
for e < 0, while a non-trivial fixed point appears for 
e > 0. As in Wilson's method, it is easy to see that 
for e < 1/12 all the higher-degree terms in S\ff\ scale 
to zero in the IR limit. Thus, the IR-dominant part 
of the action reads 



S{Q\ 



^ A,,Tr[Q,Q,]+|^Tr[Qf], (42) 



i=0 



where and m cx T ~ T^. The bare propagator A^j 
actually depends on i, j through the difference — 
X(j), where the function X(i) is defined as follows: 
given i G [0, 2*^ — 1] and its expression in base 2 



I = J2 aj2^ Ai) = ^k~i~jV. (43) 
i=o i=o 

Hence, the quadratic term of Eq. (|42p is not invari- 
ant under spatial translations. This would make any 
explicit computation of the loop-integrals, and so of 
the critical exponents, extremely difficult to perform. 
This problem can be overcome by a re-labelling of the 
sites of the lattice [l3, [3 



After relabelling one obtains that Ay depends 
on i,j just through the difference i — j, thus S[Q] 
is traslationally invariant, and the ordinary Fourier 
transform techniques [H, [l^] can be employed. In 
particular, the Fourier representation of the propa- 
gator is 



1 



2"-l 



p=0 



~2mp{i — j) 



2k 



(bir-^+m) 



(44) 

where \p\2 is the di-adic norm of p |17| . and the 
mass m has dimension [m] = 2(7 — 1. 



An interesting feature of the action (|42p is the 
fact the propagator A in Eq. ([44|) depends on the 
momentum p through its diadic norm \p\2- If we 
look at the original derivation of the recursion RG 
equation for the Ising model in finite dimension (in 
particular to the Poliakoff derivation Q), we find 
that the basic approximation was to introduce an 
ultrametric structure in momentum space: the mo- 
mentum space is divided in shells an the sum of two 
momenta of a given shell cannot give a momentum 
of a higher momentum scale cell. This has a nice 
similarity with the metric properties of the diadic 
norm, where if m, p2 are two integers, their diadic 
norms satisfies [o] bi - P2I2 < niax(|pi|2, |P2|2)- 



The field theory defined by Eq. (1421) reproduces 
the Tr[(5'^] interaction term of the well-know effec- 
tive actions describing the spin-glass transition in 
short-range [l^ and long-range [gI, H^I spin-glasses. 
Notwithstanding this similarity, the novelty of the 
HEA is that a high-order e-expansion can be quiet 
easily automatized within Wilson's method, by 
means of a symbolic manipulation program solving 
the simple RG equation ([5]) to high orders in e. This 
is not true for such short and long-range 0, [l^ [lO] 
models, where the only approach to compute the 
exponents is the field-theoretical one. Indeed, 
nobody ever managed to automatize at high orders 
a computation of the critical exponents within the 
field-theoretical minimal subtraction scheme, either 
for the simplest case of the Ising model, because 
such an automatization to high orders is not an easy 
task [U. 

The additional order of difficulty in computing v 
with the field-theoretical approach for the HEA will 
be clear in the following, as we sketch out its main 
steps. 



The field theory defined by Eq. ((42|l can be now 
analysed within the loop expansion framework. The 
renormalized mass and coupling constant are defined 



We define the one-particle-irreduciblc (IPI) 
renormalized correlation functions 

TO = m,, + Sm, (45) 
g = mf^grZg. (46) 



r^™''^(ai6iii • • •a,„6mim;ji • • -jii^jTO^s^) = 2^^ {Qi^ aibi ■ ■ ■ Q 



Tr[Ql].--Tr[Ql]) 



ipi ' 



r 



in terms of the renormalized parameters mr,gr- 
Since this model has long-range interactions, the 
field Qab is not renormalized, and Zq = 1. 
Hence, all we need to compute are the 
renormalization constants Zg,Z2 and Sm. This can 
be done by computing the IR-divergent parts of 
-p(3,o) p(2,i) ^jj^g minimal subtraction scheme 

0. In other words, one takes the IR-limit — > 0, 
and systematically removes the resulting e-singular 
parts of the correlations functions, by absorbing 
them into the renormalization constants Zg^Z^. 
The Feynmann diagrams contributing to r,'"^''^'' , rf^'^-* 
are shown in Figs. [31 HI and their singular parts are 
in the form of 1/e, 1/e^-polcs. 



the standard result [T3| that the number of integers 
p G [0, 1^ — 1] such that |p|2 = 2^-', i. e. the volume 
of the diadic sheh, is given by 2^-'+''"^, Eq. (|i5)) 
becomes 



1 



fc-i 

^1 = y 2-^-1 3 

[to^ + 5to + 2-J(2'--i)] 

^ V2-J-^ - 

^ [TO^ + 5m + 2-J(2'^-i)]^' 



(49) 



Here wc show by a simple example how the e- 
divergent part of such diagrams can be computed. 
Let's consider the one-loop expansion of F^'^'*'''. This 
is obtained by picking up the Tr[Q^]-term in the 
renormalized IPI generating functional 



2 
3! 



(47) 



n-2 



The loop-integral 



^1 = 



^ 2"-! 
2k X! 



p=0 



(48) 



is represented by the first diagram in Fig. HI 
Eq. (|48|) has a well-defined limit ioi k ^ oo. In- 
deed, thanks to the translational invariance of the 
theory, the argument of the sum in the right side of 
(|48p depends on p just through its diadic norm. It 
follows that the sum J'l can be transformed into a 
sum over all the possible values of |p|2. Indeed, using 



where in the second line of Eq. p9| the fc — !■ oo 
limit has been taken, since the series in the first line 
is convergent. By using the fact that 5m — O(g^), 
we can rewrite ([^^ as 



j^i - y2--'"-i 

f^o K + 2-^(2.-1)] 

+0(,9,2). 



(50) 



It easy to see that J^i is IR-divcrgcnt for -> 0. 
Indeed, in the limit — > the sum over j in Eq. 
(|50p is dominated by the terms in the IR region 

= \p\2 0. The js corresponding to this re- 
gion are given by the relation mr « 2"-^, and go to 
infinity as — >■ 0, yielding a divergent sum in J^i. 
In the IR region, the sum in the right side of Eq. 
(|50|) can be approximated by an integral, since the 
integrand function is almost constant in the interval 

+ 1] for large j. Setting q = 2^^ , for to^ — )• we 
have — q log 2 dj = dq, and 
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dq 



2 log 2 JO 

6e 
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rrir 

2 log 2 

rrir 

2 log 2 
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(1- 
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g2a-l)3 



(51) 



(1 



The integral in the right side of the last line in Eq. 
((5T|) is convergent for e > 0, and diverges as e — ?> 0. 
Its e-divergent part can be easily evaluated 



^1 



4 log 2 
xT II ' 



■r I 3 

1 



1 



6e 
2cr-l 



1 - 2cr 
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12elog2 



2cr. 

o,(i; 



(52) 



0(.9.'), 



where F is the Euler-Gamma function and Oe(l) de- 
notes terms that stay finite as e ^ 0. As we will 
show in the following, these finite terms will give a 
contribution to the renormalization constants at two 
loops. 

By plugging Eq. ([5^ into Eq. ([TT]) . one can com- 
pute the (^^-coefficient of Zg by imposing that the 
e-singular part of J^i is cancelled by Zg. For n ~ 
we have 



Zg = l 



1 



48elog2 



9r 



0{gt). 



(53) 



By repeating the same computation for Fr[Q, i^], 
and imposing that the X]i=o ^ Qf -term is finite, 

(2 1).. 

i. e. that FJ. is finite, we obtain 



Zo 



1 



24e log 2 



Such procedure has been pushed to two loops by an 
explicit calculation. Even if the evaluation of the 
e-divergent part of the two-loop diagrams is more 
involved, the techniques and underlying ideas are ex- 
actly the same as those used to compute the one-loop 
diagram J'l. In Figs. [31 [H we show the Feynman 
diagrams contributing to the finiteness conditions of 
Fr^'^'' and of Fr^'°'' respectively. The diagrams in Fig. 
[3]evaluated at zero external momenta will be denoted 
by J^i, . . . , J^6, while those in Fig. [4]by J'-j, Jxq. 
It is easy to show that the equalities 



^1 = 


^7, 


^2 = 


^6 = 


^3 = 


^10, 


^4 - 


^5 = 



hold, and so that all we need to compute the renor- 
malization constants are ..^i , ^2 7 -^3 , -^4 • -^i is given 
by Eq. (|48l) . while the other loop integrals are 
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In the limit m,. — > 0, J^2, ^3, ^4 are given by 
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The finiteness of vi^'^'' is imposed by making fi- 
nite the X]i=o ^ Tr[Qf]-term in the IPI generating 

(2 1) 

functional rr[Q]. The finiteness of is im- 

posed by making finite the X]?=o ^ ^''^iQf -term in the 



IPI the generating functional with Tr[Q2]_ijisertions, 
rr[Q, i^]- The two-loop expansion of rr[Q] and of 
r^Q,^] read 



r,.[Q,o] 



2*=-! 

.9. ^Tr[Q3], 

i=0 



3e 
2cr-l 
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(54) 



^ .g,m. [6(" - 2)^8 + 6(71 - 2) V- + 2(n(,i - 1) - 4 - (n - 2)2)^io] + 0{g',) \ + 



3!2 



i=0 



16 



+=^3] + 0(5')^ + 



r 



where the ■ • • in Eq. (15^ stands for terms that are 
not cubic in Qi , while the ■ • • in Eq. (j55[) for terms 
that arc not quadratic in Qi and linear in Ki. The 
renormalization constants Zg, Zi are calculated by 
imposing that the rcnormalized correlation functions 
p(3,o)^ p(2,i) finite, i. c. that the term in curly 



brackets in Eq. ([M]) and that in Eq. (f55|) have no 
singularities [7| in e. At this purpose, we observe 
that the finite part of the integral J'\ contributes 
to the renormalization constants at two loops. For 
example, let us consider the second addend in curly 
brackets in the right side of Eq. By using Eq. 

(|52|) and the one- loop result ([53| for Zg, it is easy 



to see that this term produces an e-divergent term, 
given by 



vrir 

48 



48e log 2 



(56) 



where Oe(l) is the finite part of Jx in Eq. The 
term in Eq. (|56p is of 0{g^) and singular in e. Hence, 
it contributes to the 0{g^)-term in Zg. 

After setting n = and imposing the finiteness con- 
ditions, we find 
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48e log 2 
24e log 2 
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1536e2(log2)2 
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512elog2 
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0{g'r). 



(57) 
(58) 



It is also easy to verify that Sm ^ 0{g^). 

Once the IR-safe renormalized theory has been 

constructed, the effective coupling constant .g(A) at 

the energy scale A is computed form the Callan- 

Symanzik equation in terms of the /3-function by set- 
1 

ting /i = mr"^^ 



/3(3(A)) - A 



g,m 



dgjX) 
dX ' 



(59) 



/3{gr) can be explicitly computed in terms of the 
renormalization constant Zg by applying /i 
both sides of Eq. ((iS)) 



= /i 



d_ 



(60) 



The right side of Eq. ((60|) can be then worked out ex- 
plicitly by using the two-loop result ([58|) and substi- 
tuting systematically [i[gr) to ^J^^^ ■ In this way, 

an explicit equation for P{gr) is obtained. By solving 
perturbatively this equation to 0{g^) one finds 



-3e.gr + 



9r 



8 log 2 



+ 3 



5 + 2 • 2^/3 
128 log 2 



-9r 



0{gl) 



(61) 

Setting g* = g{\ = 0), we see from Eq. (pT|) that 
the fixed point g* ~Q \s stable only for e < 0, while 
for e > a non-Gaussian fixed point g* of order 
e arises, as predicted by dimensional considerations 
and by Wilson's method. Now the IR-limit A — )- 
can be safely taken, and the scaling relations yield v 
in terms of g* and Z2 



V2[gr 



d log Z2 



r,2[g;]+2a-l 



(62) 



By plugging the two-loop result for g* and Z2 into 
Eq. ((S^ , we reproduce the result (|40p derived within 
Wilson's method. 

We observe that the analytical effort to derive the co- 
efficients of the e-expansion in this field-theoretical 
approach is much bigger than that of Wilson's 
method. Indeed, in the minimal subtraction scheme, 
additional calculations are needed to extract the 
coefficients of the poles in 1/e, 1/e^ of the Fey- 
namnn diagrams in Figs. [3l [H It follows that for 
an automatized implementation of the e-expansion 
to high orders, Wilson's method turns out to be 



much better-performing that the field theoretical 
method. Notwithstanding this, the tensorial oper- 
ations needed to compute Q-dependence of the dia- 
grams in this field theoretical approach turn out to 
be exactly the same as those needed in Section [III 
and no additional effort has been required to com- 
pute them. 



IV. CONCLUSIONS 

In a previous work (30| . we set up two pertur- 
bative approaches to compute the IR behavior 
of a strongly-frustrated non-mean field spin-glass 
system, the HEA model. The two methods yield 
the same prediction at two loops for the critical ex- 
ponent v related to the divergence of the correlation 
length. 

In this work the two-loop computation is shown 
in all its most relevant details, so that the reader 
can reproduce it. Moreover, we show the underlying 
renormalization group ideas implemented in the 
two computation methods. Amongst these, the 
existence of a characteristic length ^ diverging at 
the critical point, where the theory is invariant 
with respect to changes in its energy scale. In 
addition, we show with an explicit example that 
such a computation of the critical exponents could 
be quiet easily automatized, i. e. implemented in a 
computer program, in order to push this e-expansion 
to high orders in e, and so eventually make this 
theory physically predictive. Indeed, we give a 
graphical interpretation of the cumbersome tenso- 
rial operations needed to compute v and previously 
used in [30| . Such a graphical method makes the 
calculations much more straightforward and elegant. 
We observe that once this high-order series in e 
will be known, some resummation technique will be 
needed to make the theory predictive, because the 
series has probably a non-convergent behavior. If 
the high-orders series could be made convergent by 
means of some appropriate re-summation technique, 
this calculation would yield an analytical control 
on the critical exponents, resulting in a precise 
prediction for a non-mean field spin-glass mimicking 
a real system. 

We are glad to thank N. Sourlas, S. Franz and M. 
Mezard for interesting discussions and suggestions. 



FIG. 3: One and two-loop Feynman diagrams contribut- 
to rj.^'^'. The crosses represent Tr [Q^] insertions. 
From left to riglit, sucli diagrams computed at zero ex- 
ternal momenta are equal to ^i, J^2, =^3, ^4, ^5, ^6 re- 
spectively. 




FIG. 4: One and two-loop Feynman diagrams contribut- 
ing to Fr^'"'. From left to right, such diagrams computed 
at zero external momenta are equal to ^7 , J^g , J^g , =#10 
respectively. The last diagram is non-planar. 
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